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ABSTRACT 

We consider the dynamics of dust and gas during the clearing of protoplanetary discs. We work 
within the context of a photoevaporation/viscous model for the evolution of the gas disc, and 
use a two-fluid model to study the dynamics of dust grains as the gas disc is cleared. Small 
(< lOfim) grains remain well-coupled to the gas, but larger (~ 1mm) grains are subject to 
inward migration from large radii (~ 50AU), suggesting that the time-scale for grain growth in 
the outer disc is ~ 10'*-10^yr We describe in detail the observable appearance of discs during 
clearing, and find that pressure gradients in the gas disc result in a strong enhancement of the 
local dust-to-gas ratio in a ring near to the inner disc edge. Lastly, we consider a simple model 
of the disc-planet interaction, and suggest that observations of disc masses and accretion rates 
provide a straightforward means of discriminating between different models of disc clearing. 
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1 INTRODUCTION 

It is now well accepted that newly-formed, solar-type stars are sur- 
rounded by primordial discs of ga s and dust. Such discs can occa - 
sionally be resolved directly (e.g. iMcCaughrean & 0'DellllT99^ . 
but their presence is more commonly inferred from secondary ob- 
servational diagnostics. The most common such diagnostics are 
infrared excess, arising from the reprocessing of stell ar radia- 
by warm dust in the inner d isc disc iK enyon & Hartmannl 



' tion b y warm dust m the inner d isc disc ([K enyon & Hartmani 
1 19951 : iHaisch. Lada & Ladj I2OOII : iHartmannet al.. ,2005l) . accre 



tion signatures such as line emissi o n and ultraviolet excess 
(Hartigan, Edwards & Ghandour 1995 ; Calvet & Gullbring| ll998l : 
iGullbrin g et al. 2000 ,; MuzerpUe et al . 2000). and millimetre con- 
tinuum emission from cold dust i n the outer disc jSeckwith et al.l 
Il99d ; [Andrews & Williamj|2005l) . These different diagnostics are 
usually seen together, and few transition objects are known to ex- 
ist between the disc-bearing classical T Tauri stars (CTTs) and 
the disc-less weak-lined T Tauri stars (WTTs). Since CTTs evolve 
into WTTs, the apparent lack of transition objects implies that 
the transition between t he CTT and WTT states is very rapid 
(Skrutskieetal. 1990; K envon & Hartmannl [l99^ ; ISimon & Pratd 
ri995l ; IWolk & Waited 199^ : typical estimates suggest that the discs 
are cleared on a time-scale of order lO'yr after a lifetime of a few 
Myr. Moreover, the data show that the inner and outer discs ar e 
cleared almost simultaneously (e.g. [Andrews & Williamsl liool) . 
suggesting that the mechanism responsible for the rapid disc clear- 
ing operates over the entire radial extent of the disc, from < 0. 1 AU 
to > lOOAU. 

More recent studies have begun to provide some insight into 
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the disc clearing process, and a handful of "transitional" discs, 
with properties between the CTT and WTT states, have now been 
well-studied (e.g. GM Aur, 'Rice et al.' '2003'; 'Calvet et al.' '2004 
TW Hya, Calvet et al. 2002; Eisner, Chiang & Hillenbrand 20oi; 
CoKu TauM. IForrest et all2 004; d'Al essio et al.ii2 005a). These ob- 
jects are characterised by a lack of near-infrared excess emission, 
but "normal" disc spectral energy distributions (SEDs) at wave- 
lengths longer than ~ 10//m. However, some of these objects, such 
as GM Aur, are actively accreting and possess massive (~ O.IM0) 
discs, while others, such as CoKu Tau/4, show no evidence of ac- 
cretion and possess much l ess massive (~ IMj^p) discs. Recently, 
ISicilia-Aguilar et 'Zl j2006h identified 17 "transition objects", and 
found that roughly half showed evidence of accretion, while half 
did not. It does not, therefore, seem likely that "transition objects" 
identified in this way (i.e. on the basis of their SEDs) represent a 
homogenous class of objects. However, all such objects do seem 
consistent with some form of "inside-out" disc clearing, during 
which the disc appears to have a large central cavity or hole, at least 
in the dust disc (see also Padgett et al. 2006). The mechanism(s) 
responsible for producing such holes remain uncertain, however, 
and previo us studies have variously invoked th e presence of a 



planet (e.g. iRice el^ l2003l : lOuillen et aT]|2004l), dust evolution 



(e.g. Wilner et al. 2005) or photoevaporation (e.g. [Alexander et afl 
|2006b; 

Gpto et al. 2006) to explain the observed "holes". 

In this paper we seek to explore the behaviour of the dust 
component of the disc during the clearing phase, within the 
context of photoevaporation models of disc evolution. We also 
seek observable diagnostics that can be used to discriminate 
between different models of disc clearing. Photoevaporation 
occurs when ultraviolet radiation heats the disc surface suffi- 
ciently that it becomes unbound from the central star and flows 
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from the disc as a wind (se e the reviews by iHoUenbach et alj 
l200d iDullemond et^ |2006| . and references within). Here we 
focus on the effect of ionizing radiation from the central star, 
as the effects of external photoevaporation from nearby O-stars 
are only significant for the minority of TTs which are found 
near to the centre of massive clusters dStorzer & HollenbachI 
1 19991 : i isner & Carpentej |2006|) . It has been demonstrated that 
it is reasonable to treat the central star as having an ionizing 
luminosity of order O = 10'''photons s"', presumably arising in 
the stellar chromosphere, which is approximately constant over 
the duration of the TT phase ( [Alexander. Clarke & Pringl3 
l2005h . Most previous studies of the evolution of pho- 



toeyaporatmg 
200 ll: 
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Matsuvama. Johnstone & Hartmanri 20031: 



discs (e.g. 



Armitage. Clarke & Palld l2003l : [Alexander. Clarke & Pringlel 
2006b ') have modelled the evolution of the gas disc only, and 



simply assumed a constant dust-to-gas ratio in order to compare 
to observed data. These studies have shown that models which 
incorporate photoevaporation and viscous accretion can reproduce 
a number of observed properties, such as the rapid transition from 
the CTT to WTT state and the near-simultaneous clearing of the 
inner and outer discs. The basic premise of such models is that 
the low mass-loss rate of the wind (~ lO^'^Mgyr"') is initially 
negligible compared to the disc accretion rate. However, the disc 
accretion rate falls with time, eventually becoming comparable to 
the (constant) wind rate. At this point the wind opens a gap in the 
disc, depriving the inner disc of resupply. The inner disc drains on 
to the star on a viscous time-scale, and from this point onwards the 
wind dominates the evolution, resulting in a rapid clearing of the 
disc from the inside-out. 

iTakeuchi. Clarke & LiijfeOOl see also lTakeuchi & Lirj2005l) 
combined a viscous/photoevaporation model for the evolution of 
gas discs with a two-fluid model for the evolution of the dust, and 
studied the evolution of mm-size grains over the lifetime of the 
disc. They found that mm grains migrate inwards (due to head- 
wind drag) on fairly short time-scales, typically ~ IMyr or less. 
Consequently, in their model the mm-size grains were removed 
from TT discs prior to the dispersal of the gas disc, and they pre- 
dicted the existence of a population of dust-poor gas discs. More re- 
cently, however, Andrews & Williams ( 2005) surveyed more than 
100 TTs in the Taurus- Auriga cloud and found that essentially all 
objects (> 90%) with optically thick, accreting inner discs also 
show a strong millimetre flux. In addition, only a handful of ob- 
jects (< 10%) with no infrared excess or accretion signatures were 
detected at milli metre wavelengths. Th is tells us, in contrast with 
the predictions of ITakeuchi et al.l ( |2005|) . that mm-sized dust in the 
outer disc is dispersed at the same time as the gas disc. Moreover, 
this suggests that mm grains in the outer disc must be replenished 
throughout the lifetime of the CTT phase, as the inward migration 
time-scale of these grains is much shorter than the typical CTT life- 
time. In addition, detailed models of dust coagulation have shown 
that replenishment of grains is necessary even to sustain a popula- 
tion of small, mic ron-sized grains, over the (sev eral Myr) duration 
of the CTT phase jPuUemond & Dominikll 20051) . 

Consequently, here we make the assumption that the dust- 
to-gas ratio (for all grain sizes) is constant at the beginning of 
the clearing phase, and consider the differential motion of the gas 
and dust during the photoevaporatively-induced transition. Essen- 
tially we assu me that some replenishment pro cess(es), such as that 
suggested by iDuUemond & DominikI j2005l) . maintains a quasi- 
equilibrium in the grain population up to the point where photo- 
evaporation alters the gas dynamics significantly, and we then fol- 



low the differential evolution of the dust and gas discs through the 
transitional phase. W e use a model similar to that presented by 
ITakeuchi et al. ( 2005), but also include updated photoevaporation 
models (Font et al. 2004; Alexander et al. 2006a, b) which radically 
alter the evolution of the outer gas disc at late times. In addition we 
consider a range of grain sizes, from small /jm grains up to larger 
mm-sized particles. The structure of the paper is as follows. In Sec- 
tion|2]we describe our disc model, and present the results in Section 
[3] We discuss the implications of these results, and the limitations 
of our analysis, in Section |4] In Section [5] we compare our results 
to those from a simple model of the disc-planet interaction, and 
suggest a simple observational diagnostic that can be used to dis- 
criminate between these different models, before summarizing our 
conclusions in Section|6] 



2 MODEL 

In order to model the dynamics of dust and gas du ring the transi- 
tion ph ase, we adopt a model similar to that used bv lTakeuchi et al.l 
( l2005t) . We assume that the gaseous component of the disc evolves 
subject to viscosity and photoevaporation, and that the dust compo- 
nent is subject to diffusion and gas-drag forces. We neglect other ef- 
fects, such as grain replenishment, radiation pressure, or Poynting- 
Robertson drag; the possible consequences of these simplifications 
are discussed in Section|4] 



2.1 Gas disc 

The evolution of th e gas surface density, S„(j?, f ) is governed by the 
diff^usion equation fcynden-Bell & Pringlj|l974l ; |Pringlelll98ll) 



5Eg 
dt 



3 _a_ 

RdR 



(1) 



where t is time, R is cylindrical radius, v is the kinematic viscosity 
and the term E„ii,(](i?, t) represents the mass-loss due to photoevap- 
oration. 



2.1.1 Viscosity 

The evolution of the gas disc is determined primarily by the form 
of the kinematic viscosity v, which governs the transport of angu- 
lar momentum in the disc. The details of how angular momentum 
is transported in accretion discs are still subject to much debate 
(e.g. Gammie 1996; Balbus & Hawley jj98), but for simplicity we 
adopt an alpha-disc model dShakura & Sunvaevlll973h . The viscos- 
ity takes the form 

v{R) = aClH^ , (2) 

where a is the standard Ishakura & SunvaevI ( Il973h viscosity pa- 
rameter, n(R) = -^/gmT/W is the orbital frequency, and H(R) is the 
disc scale-height. We adopt a scale-height consistent with a "flared 
disc" model (e.g. ,Kenvon & Hartmann 1987.^ . which takes the form 



H(R) oc R'' , 



(3) 



with the power-law index q = 5/4. We normalise this relationship 
so that the disc aspect ratio H/R = 0.0333 at R = lAU. (With 
this scaling we therefore have H/R = 0.019 at R = O.IAU, and 
H/R = 0.105 at = lOOAU.) This choice of the index q im- 
plies that the disc midplane temperature scales as and re- 
sults in a viscosity which scales linearly with radiu s. Arguments in 
favour of such a viscosity law were discussed by iHartmann et al] 
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1998h. and recent observati ons also support this parametrization 
I Andrews & Williamsl |2003) . We consider two viscosity models 
within this framework: a "high-viscosity" model with a = 0.01, 
and a "low-viscosity" model with a = 0.001. 



2.1.2 Wind model 

Here we seek to model the transitional phase more accurately than 
before, and so the wind term is given by 



■ + 2d' 



(4) 



where Ejiffusc and Ejiicct represent the contributions to the wind from 
the diffuse and direct stellar radiation fields respectively. We adopt 
the diffuse wind profile from Font et ^ ( |2004|) . and the direct wind 
profile from Alexander et al. (2006aj) (with HjR = 0.05). The nu- 
merical forms of these wind profiles are defined in A ppendix IaI 
The direct profile is defined in lAlexander et alj ( l2006a l) as a power- 
law: it is only defined for R > i?;,,, and diverges to small radii. In 
order to avoid numerical problems we smooth this profile at the 
inner disc edge /?i„ by multiplying idircct(^) by the smoothing func- 
tion 



f(R)- 



1 + exp 



(5) 



where /fin is the disc scale-height at the inner edge. Both the dif- 
fuse and direct wind rates scale as the square root of the respective 
ionizing fluxes, and we model the transition by scaling these fluxes 
appropriately. The diffuse flux ©diffuse arises from recombinations 
of atomi c hydrogen in the disc a tmosphere, and was modelled in 
detail bv lHollenbach et al.l fl994h . Here we note that as the inner 
disc drains the structure of this atmosphere is modified, and the 
contribution from the atmosphere at a given radius is essentially 
removed if the disc at that radius is optically thin in the vertical 
direction to ionizing photons. We denote this radius (at which the 
disc beco mes optically thin in the vertical direction) by iJthin and, 
foUowing lHollenbach et al] h994h . we assume that the diffuse ion- 
izing flux reaching larger radii is dominated by the flux emitted at 
this radius. This flux is subject to geometric dilution, and conse- 
quently we scale the diffuse ionizing flux as 



«di 



if Rlhm < Rait 

otherwise 



(6) 



where O is the stellar ionizing flux and Ri„ is the radius of the inner 
disc edge outside the gap. (Note that Rm^ ^ Rm, by construction.) 
j?cri t is the "critical radius" at which the gap opens, which is given 
by jLilfmaijl2003l : lFont et alj2004l : lDullemond et alj200^ 



^ciit 



1.4 — ^ AU, 



(7) 



an d is approximately one-fi fth of the "gravitational radius" defined 
bv lHollenbach et aU ( Il994l ). If the surface density in the inner disc 
is sufficiently small that the disc is optically thin to ionizing photons 
i^thin is located where 



2g(/?(hin) 



(8) 



Here (Tu f^y = 6. 3 x 10~'^cm ^ is the absorption cross section for 
ionizing photons ( |Cox||200^ . and is the mass of a hydrogen 
atom. 

The direct radiation field reaching the outer disc is simply the 
stellar field attenuated by material in the inner disc. Consequently 
we evaluate Odirect as 



O exp(-T) , 



(9) 



where t = Nau ^ev is the radial optical depth to ionizing pho- 
tons due to attenuation by gas with radial column density A'. We 
neglect attenuation due to dust in the inner disc (see Section l3Tt . 
[Alexander et al.l ( [2006ah found that most of the mass-loss due to the 
direct field comes from 1-3 scale-heights above the disc midplane, 
so we evaluate the gas column A' as 



N : 



n^dR , 



(10) 



where the gas number density n^{R) is evaluated from the surface 
density as 



nAR) 



■ exp(-;r'/2) . 



yl2nH{R)mn 



(11) 



Here the exponential term accounts for evaluating the density at x 
scale-heights above the midplane. We adopt = 2.5, but note that 
the evolution is not especially sensitive to the choice of x- 

Consequently the transition from the diffuse t o direct regime is 
treated more accurately than in previous work (e.g. lAlexander et aU 
l2006b l). At early times the inner disc is extremely optically thick to 
ionizing ph otons, so the wind profile is simply the "diffuse" pro- 
file of Fo nt et alj l l2004h . Similarly at late times the inner disc is 
optically thin to ionizing p hotons at all radii , and th e wind is sim- 
ply the "direct" profile of I Alexander et all ( l2006ah . A more de- 
tailed approach could consider the progression of the ionization 
front through the (draining) inner disc, using something akin to a 
Stromgrem criterion. Evaluating the attenuation in this manner re- 
sults in less attenuation than the parametrization described above, 
and leads to a slightly earlier rise in the direct field. However, this 
form is also rather sensitive to the choice of the disc scale-height, 
and the difference between the two parametrizations is not signif- 
icant. (The direct field "turns on" slightly earlier, but the differ- 
ence is only a small fraction of the already short viscous time-scale 
of the inner disc.) Consequently we adopt the more conservative 
parametrization described in Equations 16141 II 



2.1.3 Initial conditions 

We are only interested in the dynamics during the transition from 
CTT to WTT, so we choose initial conditio ns that reflect this 
fact. Previous p hotoevaporation models (e.g. IClarke et al.l l200ll : 
[Alexander et alj l2006b) have shown that the wind starts to influ- 
ence the evolution when the disc accretion rate drops to around 1/4 
of the integrated (diffuse-field) wind rate, so we choose initial con- 
ditions to mimic this stage in the evolution of the gas disc. We adopt 
an ionizing flux of O = lO^'^photon s"' throughout. This results in 
a (diffuse-field) wind mass-loss rate of ^ 4 x lO^'^Moyr"' 
choose an initial accretion rate of Mq = 1.0 x 10""'MQyr 
defi ne the initial s urface density profile to be that of a steady disc 
(e.g. |PringIelll98ll) : 



, so we 
We 



i:ar) 



Mo 




(12) 



Here the second term is due to the presence of a (zero-torque) inner 
boundary at /?i„, and the exponential cutoff beyond radius R^ de- 
fines the size of the disc (and prevents divergence of the disc mass 
to large radii). We adopt R^ = 500 AU throughout. 
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2.2 Dust disc 

In order to model the dynamics of tlie dust component of the disc, 
we follow the approac h adopted bv lTakeuchi & Lij l l2005l see also 
iTakeuchi & Lirj|2002 l). We make the simplifying assumption that 
the grains are not replenished over the short time-scales in which 
we are interested, so the only way that dust mass can be "lost" is by 
accretion on to the star. The evolution of the dust surface density 
l-iiR, t) is governed by the continuity equation 



dLi 1 d 
at R oR 



(13) 



Here Fjjf is the diffusive mass flux, and is the radial dust veloc- 
ity induced by gas-drag. Our initial conditions result in an outward 
pressure gradient in the gas disc, which causes the gas at a given 
radius to orbit at a slightl y sub-Keplerian velocit y. This gives rise 
to a headwind drag force ( Weidens chillingiri97 7^. which removes 
angular momentum from the dust and results in inward radial mi- 
gr ation of the grain s . The form of this velocity was studied in depth 
bv lTakeuchi & Li j j200 j) . and we adopt their prescription for : 



(14) 



Here mk = (GM,/R)^'^ is the Keplerian orbital velocity, and the gas 
radial velocity «„ is given by 

3 d 



R^/^1.., clR 



(15) 



Ts is the non-dimensional stopping time due to gas-drag on grains 
at the disc midplane (i.e. the stopping time divided by the orbital 
time), and is given by 



(16) 



for grains of radius s and density pj. We fix the density to have the 
fiducial value pd = Ig cm"^, and explore various values of s. rj is 
defined as as the ratio of the pressure gradient to the gravity, and is 
given b}^ 



//\2/dlogEg 



+ (? - 3) 



(17) 



.Rl \ dlogR 

The diffusive term F^jf was derived by IClarke & Pringld 

( Il987l) . for the case where Sa/Sg <K 1, and in this case is given 
by 



d /S, 



(18) 



The relationship between the dust diffusion coefficient D and the 
gas difi^usion coefficient (viscosity) v is not fully understood, with 
hydrodyna mics models suggesting ratios of v/D t hat range from 
^ 0.5-10 JCar ballido. Stone & Pringle 2005; Joha n^en & KlahJ 
I2OO5I ; [Turner et alj|200a : IFromang & Papaloiz ou 2006). For sim- 
plicity we assume that the dust difi^uses in a comparable manner 
to the gas, and therefore adopt J5 = v, As long as the concentration 
Sd/Sg ^ 1, or as long as the advective term dominates the dust evo- 
lution, this form for the diffusive term is valid. However, our model 
breaks down at late times, as the dust-to-gas ratio at the inner edge 
of the gas disc can become large. Given the other uncertainties in 
the model, we make no attempt to correct for this and simply stop 
the calculations at this point. This is addressed Section |4!2l 



' Note that our definition of the power-law index cj dift'ers from that adopted 
bv lTakeuchi & Lij i2005h ; Equation ll7l has been adjusted accordingly. 



The ability of the photoevaporative wind to carry dust grains 
depends on whether gas drag from the wind can ov ercome gravity 
and lift the grains away from the disc. As shown bv lTakeuchi et al.l 
(2005), for spherical grains in a laminar flow the drag force from 
the wind is approximately 



(19) 



where is the mass of a grain, p„ is the density of the gas at 
the base of the flow, and the flow is launched at the sound speed 
Cs = 10km s '. The wind densi ty peaks close to the in ner edge of 
the disc in the direct field case ( [Alexander et al.ll2006al) . and if we 
compare F„ to the weight of the grains, F„ = GMjiiilR^, we see 
that that F„ < Fg as long as ^ > l//m. Moreover, we note that dust 
grains are expected to sediment towards the disc midplane as the 
disc evolves, but the photoevaporative wind typically flows from 
several scale-heights above the disc midplane. Thus only small, 
sub-;um sized grains can be carried by the wind, and only then if 
something (such as turbulence) lifts them to the upper regions of 
the disc. Consequently we neglect dust mass-loss due to the photo- 
evaporative wind. 

We define the initial conditions such that Xd/£g = 0.01 at all 
radii. Our model considers only a single grain size s, so formally 
this implies that a single grain size has a dust mass that is 1% of the 
gas mass. However, we note that the evolution of the dust disc in 
independent of the absolute value of Sd> as long as the dust-to-gas 
ratio remains small. 



2.3 Numerical method 

We solve the diffusion equation for the gas surface density (Equa- 
tionfD numerically using a stand ard first-order explicit scheme (e.g. 
IPringle. Verbunt & Waddll986l) . on a grid of points equispaced in 

We adopt zero-torque boundary conditions (i.e. we set = 
in the boundary cells). 

The equation for the evolution of the dust is evaluated simul- 
taneously on the same spatial grid. We use a staggered grid, with 
scalar quantities evaluated at at zone centres and vectors (veloci- 
ties) at zone faces. The integration is explicit, using a first-order 
method for the diffusive term and a second-order (van Leer) method 
for the advective term. Again we adopt "outflow" boundary condi- 
tions by setting Ed = in the boundary cells. However, the inner 
boundary condition introduces an artificial inward pressure gradi- 
ent in the gas, resulting in artificial outward advection of dust near 
the boundary, so we neglect the pressure gradient term in Equation 
|17|in the region close to the inner boundary. A similar, but phys- 
ically real, inward pressure gradient is introduced by gaps in the 
gas disc which form as a result of photoevaporation. However, the 
numerical form of the gap can result in rather large radial dust ve- 
locities, so in such cases we limit the magnitude of the (outward) 
dust velocity to be less than or equal to that of the (inward) gas 
velocity. 

The inclusion of the photoevaporative wind at a stage in the 
calculation where the dust surface density is not negligible in- 
troduces an ad ditional numerica l comp lica tion compared with th e 
calculations of ITakeuchi & LiiJ ( l2005h or ITakeuchi et al.l ( |2005|) . 
namely that the concentration term 2d/Sg diverges as the gas sur- 
face density becomes small (at the gap in the gas disc). In regions 
with no gas, the dust remains at a constant radius on a Keplerian 
orbit, so in cells where Eg = we set ffL^ldt = 0. However the 
divergence of the concentration term also causes a problem in cells 
adjacent to those where Eg = 0, as our evaluation of ffL^ldt uses 
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"three-cell" derivatives for both the advective and diffusive terms. 
In this case we revert to a first-order (donor cell) derivative for the 
advective term, and evaluate the diffusive term one cell further from 
the "gap". 

This procedure successfully removes the troublesome diver- 
gent cells from the numerical calculation, but at the cost of reduced 
accuracy close to the edges of the gas disc. The treatment of the 
diffusive term in particular causes problems, as the surface density 
of the gas rises sharply near the gap, and the accuracy is limited 
by the grid resolution. Essentially we evaluate the diffusive term 
one cell-width away from where we should, so lowering the grid 
resolution reduces the accuracy still further. The consequence of 
this inaccuracy is that dust mass is not exactly conserved close to 
the gap in the gas disc. In order to minimize the effect this has on 
the overall accuracy of the calculation we take a "brute force" ap- 
proach and simply use a large number of grid cells. We adopt a grid 
spacing of AR''^ = 0.0125AU''^ and model the disc over the radial 
range [0. 1 AU,2000AU], so consequently we use 7108 grid cells. At 
this resolution, dust mass is conserved to better than ±5% in all our 
models. 



3 RESULTS 
3.1 Gas disc 

We first describe the evolution of the gas disc alone. The evo- 
lution does not differ s ignificantly from the models presented in 
I Alexander etai] ( l2006bh . and we refer the reader to this work for a 
detailed description of the gas disc evolution. The evolution of the 
gas surface density Eg in the high- and low- viscosity models, can 
be seen in Figs|2]&[3]respectively. Our models begin (i.e. they have 
r = 0) at the point where the wind begins to influence the evolution, 
so the wind opens a gap in the disc rather rapidly, on the viscous 
time-scale of the inner disc. Once the surface density in the inner 
disc drops below Eg ^ 10^'g cm"-, the inner disc becomes optically 
thin to ionizing photons and the direct field begins to dominate the 
evolution, rapidly clearing the outer disc. 

By comparing the two viscosity models, we see that the abso- 
lute clearing time-scale varies almost exactly as 1 la. In the diffuse 
wind regime the evolution is driven by the viscosity, so the time- 
scale to open the gap scales exactly as l/or. Once the direct field 
takes over, however, the wind dominates the evolution (rather than 
the viscosity, see discussion in lAlexander et al. ''2006b) and we see 
a small departure from this l/o- scaling. This occurs because the 
gap is opened at constant accretion rate, so the low-viscosity disc 
is more massive and the wind takes correspondingly longer to dis- 
perse it. 

In the high-viscosity model, with a = 0.01, the initial mass of 
the gas disc is 4.6 x IO^^Mq, or around 0.48 Jupiter masses. The 
effect of the wind on the disc is apparent almost immediately, and 
a well-defined gap in the disc appears at f = 3.3 x 10"* yr. Cut-off 
from resupply, the inner disc then drains on its viscous time-scale. 
In the low-viscosity case (a = 0.001), the initial disc mass is larger 
by a factor of 10 (4.6 x lO^^Mg, or around 4.8 Jupiter masses), 
and the evolution slower by a factor of 10 (i.e. the gap opens at 
f = 3.3 X lO'yr). As the inner disc drains it eventually becomes 
optically thin to ionizing photons, and the direct field begins to clear 
the outer disc: this occurs at f = 4. 1 x lO^'yr for a = 0.01, and 
4.3x lO^yr for a = 0.001. We stop the calculations at / = 5.5x lO'^yr 
(5.5 X lO^yr for the low-viscosity model), at which point the inner 
edge of the disc has reached approximately 60.5AU (47.5AU). 



o 




O : 



10"^ 0.01 0.1 1 10 100 1000 

s / mm 

Figure 1. Orbital decay times, defined as f = -R/u, from R = lOAU in 
the initial, steady, disc, plotted as a function of grain size .s. The heavy 
hnes show the decay time due to advection, while the light lines show the 
(constant) viscous time-scale for diffusion of the gas. The solid lines show 
the results for the high- viscosity model; the dashed lines the low-viscosity 
model. 

The more accurate treatment of the diffuse-direct transition 
(described in Section [2.I.2t re sults in a smooth e r transi tion between 
the two regimes than seen by [Alexander et al ] l l2006bl) . but the dy- 
namics are not significantly modified. As we show in the following 
Section, the dust in the inner disc is in fact removed more rapidly 
than the gas, so the inclusion of dust opacity has no effect on the 
results. As mentioned above, a more accurate treatment still could 
consider the evolution of the ionization front through the draining 
inner disc, using something akin to a Stromgren criterion, but this 
should not result in significantly different behavioui^. In this case 
the inner disc will become fully ionized when the surface density 
drops below ^ 10""'g cm"- (the exact value depends rather sensi- 
tively on the disc thickness, as the recombination rate scales with 
the square of the gas density), so we expect that the direct field will 
become significant once the total gas mass in the inner disc falls 
below ~ IO^'Mq. This results in a the direct field becoming sig- 
nificant slightly earlier than in our model, but the difference is not 
significant for the evolution of the disc. Moreover, given the other 
uncertainties in this process, such as the behaviour of angular mo- 
mentum transport in the low-surface density inner disc, we consider 
this treatment to be sufficiently accurate. 

3.2 Dust disc 

In order to understand the evolution of the dust disc, it is instructive 
to re-cast Equation [13] in terms of a diffusive radial dust velocity 
Uiif, so that 

^ + [RI,i(u,,! + u,)] = . (20) 
ot R oR 

The advective velocity is given in Equation [14] and depends on 
both the grain size and gas density, while we can express the diffu- 
sive velocity as 

^ This is done by solving the recombination balance equation <I) = 
J asn^dV to find the gas number density n. where dV is the volume ele- 
ment and ffB is the appropriate (Case B) recombination coefficient. 
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Mdif ■ 



£> rflog(£j/£g) 
R dlogR 



(21) 



The diffusive dust velocity depends on the viscosity, and also the 
dust concentration gradient, ii^n is therefore comparable to the dif- 
fusive gas velocity Hg, but is usually slightly larger because of the 
influence of the concentration gradient term. The diffusive com- 
ponent acts to smooth out the concentration gradient, while the 
advective term always acts radially inward (until the gap in the 
gas disc opens). We see, therefore, that the evolution of the dust 
grains is determined by competition between these two compo- 
nents; whichever acts more rapidly dominates the evolution. 

Fig[T] shows how the orbital decay time-scale (at R = lOAU) 
due to advection, t = -Rju^, varies as a function of grain size s, 
in the initial disc. Also shown is the (constant) orbital decay time- 
scale for the gas disc (essentially the viscous scaling time ty). As 
we have adopted D = v, the diffusive time-scale of the dust disc is 
simply given by 



fdif - 



rflog(£j/£g) 



(22) 



We see, therefore, that the dust diffuses on a time-scale that is typi- 
cally a factor of a few shorter than the viscous time-scale. However, 
if the concentration gradient becomes large, then the dust diffuses 
much faster. In our models the concentration gradient term is al- 
ways greater than unity, but rarely becomes larger than ~ 5 ex- 
cept close to the gap in the disc. We therefore see from Fig[T]that 
(prior to the opening of the gap) the evolution of small grains, with 
s < lOyum, is dominated by the diffusion, while the evolution of the 
larger grains is dominated by advection. Thus we expect the con- 
centration of small grains to remain approximately uniform up to 
the point where the wind opens the gap in the disc. Larger grains, 
by contrast, migrate inwards more rapidly than the gas disc evolves, 
so in this case we expect significant evolution of the dust inside 10- 
20AU before the gap is opened. 

Fig[T]also highlights a key difference between dust evolution 
in the high- and low-viscosity models. In the high-viscosity model 
the gap in the gas disc opens after / = 3.3 x lO'^yr, but this takes 
10 times longer in the low-viscosity model. The difference in the 
decay times for large (~ 1mm) grains, however, differs only by a 
smaller factor (around 1-5). Consequently, we expect to see a more 
pronounced inward migration of such grains in the low-viscosity 
model. 

Fig[2] shows the evolution of the four high- viscosity modelfl 
We see some inward migration of the larger grains at early times 
(before the gap opens), but we do not see a significant depletion 
of grains beyond ~ lOAU during this phase. Once the gap opens, 
however, we see a very different behaviour. The gas surface den- 
sity increases sharply outside the gap (i.e. at the inner edge of the 
outer disc), resulting in a strong inward pressure gradient. The gas 
close to the gap edge is therefore super- Keplerian, and we see rapid 
outward migration of the grains. Consequently, once the gap has 
opened the inner dust disc is cut off from resupply, and grains of 
all sizes drain rapidly onto the central star. As the direct field pho- 
toevaporates the outer gas disc, the inner edge of the disc moves 



' Note that we assume an initial dust-to-gas ratio of 0.01 in each model 
(i.e. for each grain size), rather than a total dust-to-gas ratio of 0.01. The 
absolute values of coiTespond to the surface densities we would see if 
all the dust mass was contained in grains of a particulai' size: we plot the 
different models together in order to ilfustrate the differential migration of 
different grain sizes. 



steadily outwards. During this phase the steep pressure gradient at 
the inner edge continues to drive rapid outward migration of the 
grains, with the consequence that grains of all sizes are "swept up" 
by the moving edge of the gas disc. As this process continues, more 
and more dust is accumulated into the region close to the edge of 
the gas disc, and eventually the dust-to-gas ration becomes large. 
At this point our method, treating the dust as a trace contaminant, 
is no longer valid, so we stop our calculation at this point. The "fi- 
nal" configuration is a roughly constant dust-to-gas ration at large 
radii for all grain sizes, with a concentrated dust ring at the inner 
disc edge. 

The evolution of the low-viscosity models (as seen in Fig[3]l, 
however, is rather different. In this case the larger, advectively- 
dominated grains (lOOyum and 1mm) undergo much more signif- 
icant radial migration before the gap opens. Consequently, at the 
point where the gap opens, the outer disc is significantly depleted 
in large grains (out to radii > lOOAU), by an order of magnitude or 
more compared to the high-viscosity case. Once the gap opens we 
again see rapid draining of dust from the inner disc, and again the 
inner edge of the gas disc sweeps up the dust as it moves outwards. 
However the "final" state of the evolution shows a marked deficit 
of large grains, by comparison to the high- viscosity model. We also 
note that the smaller yum-sized grains undergo significant inward 
migration from the outermost radii (> 500AU). This is due to the 
steepening pressure gradient beyond the "disc radius" R^ = 500AU, 
however, and is not physically significant (as it depends only on the 
choice of Ra). 



4 DISCUSSION 

4.1 Observational appearance of clearing discs 

A significant result from these models is that we are now able 
to make detailed predictions regarding the observable properties 
of so-called "transitional discs", which are "caught in the act" 
of clearing their discs. Our models show that the dust-to-gas ra- 
tios for small grains are approximately constant (away from the 
edge of the gas disc), so our predictions regarding near- and mid- 
infrared emission agree well with previous gas-only models (e.g. 
[Alexander et al. 2006b). The gap is opened at ^ 1 AU, and once the 
inner disc is drained the inner disc edge moves outward in radius. 
For the viscosity law adopted here (v oc R) the time required to 
double the hole radius scales as i?''^, but a steeper viscosity law 
will result in a flatter distribution of hole sizes (see Section 2 of 
[Alexander etai][2006b) . Given the short duration of the clearing 
phase relative to the disc lifetime and the corresponding paucity of 
transition discs, it is unlikely that this weak bias towards larger hole 
sizes will be observationally significant. Consequently, we predict 
that the distribution of inner hole sizes should be approximately 
uniform, from a minimum of ^ lAU outwards (for a IM© star). 

As most of the disc mass resides at large radius, the disc 
mass remains approximately constant during clearing (as seen in 
Alexander et al. 2006b). The exact value depends on the model pa- 
rameters (viscosity, stellar ionizing flux and outer disc radius), but 
disc masses of a few Jupiter masses are typical. We note, however, 
that disc masses are typically measured from dust emission at mil- 
limetre wavelengths, so these observations are subject to the uncer- 
tainties discussed in Section |431 

The constraints of the photoevaporation model require that the 
accretion rate on to the star be essentially zero during the inner 
hole phase, as once the inner disc is drained no accretion can per- 
sist. During the inner disc draining, the accretion rate falls from 
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Figure 2. Snapshots of the evolution of the gas and dust discs for the high-viscosity (a = 0.01) models: snapshots are plotted at / = 0, 2.0, 3.3, 3.6, 4.0 & 
5.0 X IC^yr. The gas surface density is plotted in black, with the dust surface density in red: different line styles correspond to different grain sizes. Note that 
these figures show the results of 4 independent models (with different grain size s), and are plotted together merely for comparison. 
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Figure 3. As Fig|2] but for the low-viscosity (a = 0.001) case. Here the snapshots are plotted at f = 0, 2.0, 3.3, 3.6, 4.0 & 5.0 X lO'yr. 



lO^'^Mgyr"' to < lO^'^Mgyr"' on a very short, viscous time- 
scale. Consequently we predict that few, if any, objects should have 
measurable accretion rates below the photoevaporative wind rate 
of ^ lO^'^Moyr"', and note that the observational detection of this 



accretion rate threshold would provide strong support for the pho- 
toevaporation scenario. 

Moreover, once the gap appears, the model predicts (and in- 
deed requires) that the inner hole region, inside ^ lAU, be es- 
sentially devoid of gas. As noted in Section ISTI in order for the 
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inner hole to be optically thin to ionizing photons the total (neu- 
tral) gas mass inside lAU must be < lO^'Mg. Recent observa- 
tions with the Spitzer Space Telescope have placed strong upper 
limits on the gas mass i n systems with cleared or clearing discs. 
iHollenbach et alj jlOOSt) looked at the nearby star HD105, and 
placed upper limits of < O.O lMjm, on the gas mass inside lAU. 
More recentlv. lPascucci et al. observed 12 nearby stars with 

optically thin dust discs and concluded that there was < O.OOOSMjup 
(< 8x 10^'Mq) of gas in the region from 0.3-1 AU. These upper lim- 
its are rather close to the threshold value above, suggesting that the 
inner disc does indeed become optically thin to ionizing photons 
at late times, and adding further support to the photoevaporation 
model for disc clearing. 

It is less clear, however, how the appearance of inner 
holes should scale with stellar mass. In recent years, obser- 
vations of disc properties, such as accretion rates and disc 
masses, have expanded to cover a much wider range in stellar 
mass than was previousl y possible (e.g. iMuzerolle et al.l 1 20051 : 
IScholz. Javawardhana & Wood2006') , and a number of trends with 
stellar mass are now known. It has been suggested that the frac- 
tion of "transition discs" (identified from near- and mid-infrared 
obse rvations of SE Ds) incre ases with decreas ing stellar mass: 
both lMcCabe et al.l (2006) and Lada et all ( l2006h found that "SED 
holes" were more prevalent among later type stars. The degree to 
which these results are influenced by selection biases remains un- 
clear, as these data are only sensitive to "holes" at < 20yum, but we 
note in passing that such an elfect occurs naturally in our models 
if the viscou s time-scale increases with de creasing stellar mass (as 
suggested bv lAlexander & Armitagell200^ . 

4.2 Dust at the inner disc edge 

An interesting result of our model is the behaviour of the dust close 
to the inner edge of the gas disc. As outlined above, once the gap 
opens the dust in the inner disc is rapidly drained onto the star. 
Outside the gap, however, the inward pressure gradient results in 
outward migration of the dust grains, and inward accretion of dust 
is halted at, or close to, the peak in the surface density profile. Con- 
sequently the inner disc cannot be resupplied with dust, resulting in 
a very "clean" inner hole. 

Once the inner gas disc is drained, the disc edge begins to 
move outward. The steep pressure gradient at the edge causes the 
grains to be "swept up" as the gas edge moves outwards, resulting 
in a rather concentrated, dust-rich, ring immediately outside the in- 
ner edge of the gas disc. As seen from Figs [2] & [3] the behaviour 
of this dust ring is largely independent of grain size over the range 
we consider. The exact form of this ring in our models depends 
rather strongly on the wind profile at the inner disc edge (defined 
by the smoothing term in Equation O, but such a ring occurs for 
any choice of parameters. Moreover, similar structures, caused by 
the same physical process (a reversal in the direction of the pressure 
gradient close to the disc edge), have been found by other authors 
(e.g. Paardekooper & Mellema 2004, 2006), so although we are un- 
able to make detailed predictions we are satisfied that the formation 
of the dust-rich ring is a real physical process. 

As the disc edge sweeps outwards the dust-to-gas ratio close to 
the edge rises, and at late times can exceed unity. At this point our 
two-fluid treatment of the disc breaks down (so we stop the mod- 
els at this point), and the evolution of the system becomes rather 
uncertain. Beyond this point the dust dynamics will be largely in- 
dependent of th e gas disc, and it may well be that such a stru cture 
is unstable (e.g. lOoldreich & Wardlll973l : lGaraud & Liij2004h . Al- 



ternatively, the enhanced dust concentration may result in enhanced 
grain collision rates, and may facilitate rapid grain growth through 
coagulation. Again, the details are somewhat model-dependent but, 
given that gas is "lost" to the wind while dust is merely redistributed 
radially, it is straightforward to estimate at what radius the dust-to- 
gas ratio of the ring will reach these large values: this typically oc- 
curs at ~ 50-lOOAU. Other studies have similarly found that local 
pressure maxima can result in sig nificantly enhanced dust concen- 
trations ( iKlahr & L irJl200ll. [2005 : Rice et al. 2006a; Besla & Wi] 
l2006l) . and speculated as to the possible consequences in terms 
of planetesimal and debris disc formation. Detailed study of these 
consequences are not possible here, but we note that the formation 
of the dust-rich ring during the gas clearing phase can in principle 
lead to rapid growth of planetesimals, and may result in the forma- 
tion of a debris disc at a radius of ~50AU. 

Were we to extend our analysis to larger grain sizes, however, 
it seems likely that grains larger than some critical size would be so 
weakly coupled to the gas that they would not be swept up in this 
manner. Inspection of Fig[T] suggests that this critical size is likely 
in the range 10- 100cm, but this is difficult to quantify without mak- 
ing detailed calculations. "Filtration" of the dust in this manner is 
analogous to the process considered by iRice et al. (2006b). They 
consider a different situation, namely the flow of gas and dust past 
a planet-induced gap at a fixed radius, but the underlying physical 
processes are the s ame. The key diffe rence between the the situa- 
tion considered bv lRiceetal] ( l2006bl) and that considered here is 
that accretion of gas persists across a planet-induced gap, whereas 
a photoevaporative wind prevents the inward accretion of gas at the 
disc edge. Consequently, in the case of a planet-induced gap small 
grains (which are strongly coupled to the gas) are carried across 
the gap by the gas, while in our models these grains pile up close 
to the disc edge. Consequently, it seems that dust properties can 
be used to discriminate between different mechanisms for creating 
"inner holes": planet-induced gaps tend to populate the inner hole 
with small, micron-sized grains while filtering out larger particles, 
while photoevaporation removes all but very large grains (cm-size 
or larger) from the inner disc. 



4.3 Grain growth and replenishment 

A significant simplification in our model is the fact that we ne- 
glect grain collisions. Collisions between grains provide a means 
of altering the size distribution of grai ns, both through coagulation 
and destructive collisions (e.g. lDominik et al. 2006) . Models have 
shown that the time-scales for grain growth and destruction can be 
significantly shorter than disc lifetimes, and have also shown that 
grains of all sizes must be subject to continuous resupply in order to 
sustain t he observed grain populations over the lifetime of th e gas 
disc (e.g. lTakeuchi et alj|200a : iDullemond & Dominikll2005l) . Re- 
cent observations have shown evidence fo r dust settling and mod- 
erate grain evolution over disc lifetimes iFurlan et alj|200p ). but 
wholesale changes in the grain population do not occur. 

Our model considers only the short, clearing phase of the evo- 
lution, but even over these short time-scales we see significant in- 
ward migration of 100pm- and mm-size grains from large radii 
(> 50AU). Observations require that this grain population be sus- 
tained even during the clearing phase, and this enables us to make 
crude estimates of the time-scale of grain growth in the outer disc. 
Our high-viscosity model shows only moderate depletion of mm 
grains in the outer disc, but the depletion is much more significant 
in the low-viscosity model. This suggests that the time-scale for the 
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Figure 4. Evolution of dust-to-gas ratio of i = 1mm grains, normalised to 
the initial value, for various values of ?rep- The dot-dashed, dotted, solid and 
dashed lines are the results for (rep = lO^yr, 10'*yr, lO^yr and lO'yr respec- 
tively. The heavy solid line shows the results for ?,-ep = 200;r/n (i.e. 100 
local orbital periods), which gives ;,ep = lO^yr at S = lOOAU. We see that 
replenishment time-scales in the range lO^-lO^ yr are required to sustain an 
approximately constant dust-to-gas ratio over Myr time-scales. 



replenishment of the mm grain population at radii of ~ 50-lOOAU 
is < lO^r. 



4.3.1 Simple replenishment model 

In order to quantify this analysis, we have considered a model 
which includes a simple replenishment term. We assume that the 
dust mass for a given grain size can be increased by growth pro- 
cesses, and add a source term to Equation [T3] The source term 
(which appears on the right-hand side of Equation 1131 ) takes the 
form 



(23) 



where f^ep is the time-scale for replenishment. We consider a model 
of the entire disc lifetime, rather than just the clearing phase, with 
parameters a = 0.01, Ra = 50AU, and Mq = 5 x lO^'^Moyr"'. 
With these parameters the accretion rate falls to lO^'^Mgyr"' after 
4.8Myr, and the gap in the gas disc begins to open at 5.2Myr. 

Fig|4] shows the evolution of the dust-to-gas ratio for models 
with s = 1mm and various values of f,.ep- We see that a replenish- 
ment time-scale of lO'yr results in an almost constant dust-to-gas 
ratio over the Myr lifetime of the disc. Longer replenishment time- 
scales result in significant depletion of dust, while significantly 
shorter replenishment time-scales result in a dust-to-gas ratio which 
increases with time. This provides an additional lower limit to the 
replenishment time-scale for mm grains in the outer disc (~ 50- 
lOOAU), and tells us that these grains must be replenished on a 
time-scale of lO'^-lO^yr, in broad agreement with models of colli- 
sional grain growth (e.g. Dullemond & Dominik 2005). Moreover, 
this provides an a posteriori justification of our choice of initial 
conditions, and suggests that replenishment processes can indeed 
sustain a quasi-equilibrium grain population over the lifetime of 
the T Tauri phase. 



4.4 Limitations 

In addition to the issues discussed above, there a number of lim- 
itations to our analysis which bear consideration. One potentially 
important physical mechanism which is neglected here is the effect 
of radiation on the dust grains, which in principle could result in ei- 
ther Poynting-Robertson (PR) drag or radiation-pressure blowout. 
However, both of these elfects are only usually significant in the ab- 
sence of a gas disc, as the gas acts both to stabilise the dust against 
radiation forces, and to shield the grains from the radiation field. In 
our models the dust is always co-located with the gas disc, so PR 
drag is negligible. Close to the inner edge of the disc radiation pres- 
sure may become significant, but the likely outcome of this process 
is the either radiative blow-out of the grains o r the formation of a 
dust-rich ring close to the edge of the gas disc taahr & Lin"200lh. 
Our models already form such structures, due to the pressure gra- 
dients in the gas disc, so it is unlikely that the inclusion of radiation 
pressure will modify our results significantly. 

In addition to these physical limitations, there are also numer- 
ical limitations to consider. Our treatment of the dust close to the 
disc edge is rather crude (see Section [23t , but it is sufficiently ac- 
curate for situations we consider. More significant, however, may 
be the form of the inner boundary condition. We adopt zero-torque 
boundary conditions throughout, but note that these can result in 
outward migration of the grains near to the boundary (see Section 
|2.3I >. In order to minimise the influence of the boundary condi- 
tion, we neglect the pressure terms in the region close to the in- 
ner boundary. The behaviour of discs close to the inner boundary 
is not well understood, but some physically-mot ivated boundary 
condi t ions, such as magnetospheric trunca tion (e.g. lGhosh & Lama 
1 19781 : iHartmann. Hewett &Calvet|[T994h . do result in conditions 
not dis-similar to those adopted here. Consequently, in such cases 
dust grains may indeed "pile up" in the region close to the inner 
boundary. We note, however, that in CTTs the inner edge of the gas 
disc usually lies inside the dust sublimation radius l ld'Alessio et al.l 
l2005bl : [Eisner et alj|2005l) , so this effect is probably not significant 
in TT discs. 

An interesting extension of this work would be to consider a 
range of grain sizes simultaneously. Such a model would enable a 
more detailed analysis of the effects of grain growth and destruction 
than are possible here. It would also allow us to consider the effects 
of factors such as the grain size distribution, and to study of the 
evolution over longer time-scales than those considered here. This 
is obviously beyond the scope of this investigation, but provides an 
interesting avenue for future work. 



5 DISCRIMINATING BETWEEN MODELS 

In light of these results, we now seek an observational diagnos- 
tic that will allow us to discriminate between the various different 
models for "inner hole" discs. We propose that the simplest such 
diagnostic lies in the relative values of the disc mass and stellar 
accretion rate. In order to compare to the results of our photoe- 
vaporation models, we now consider a simple model for a planet 
embedded in a disc. 

Several studies have shown that a relatively massive 
planet can pr oduce an observable "hole" i n a protoplane- 
tary d i sc (e.g. iRice et alj l2003l : lOuillen et all |2004| : iRice et all 
l2006bl : IVamiere et alj |2006|) . In these models the tidal inter- 
action between the planet and the disc opens a gap in the 
disc, and the tidal baiTier limits accretion across this gap. In 
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gene ral this only happens if the planet mass is > lMju„ 
(e.g.lArtvmowicz & Lubowlll996l;lTakeuchi. Mivama & Linlll996l : 



iLubow. Siebert & Artvmowicjll999f) . as less massive planets are 
not capable of opening such a gap unless the disc is unusually invis- 
cid. Moreover, the accretion rate past the planet is a strong function 
of the planet mass. 

Our "toy" model is defined as follows. We assume that the 
disc is in a steady state, and adopt the same viscosity law as used 
previously (see Section |2.1.1t . The total disc mass is therefore 
given by 



Ma 



f 



2nRI,(R)dR , 



(24) 



where Ro^ is the outer disc radius. In a steady disc, the surface 
density can be expressed i n terms of the steady disc accretion rate, 
M, and the viscosity thus ( |Pringleiri98lh 
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We can integrate Equation [24] to find 
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and we note that Q.H^/R is constant for our choice of the flaring in- 
dex, q. (A different choice of q, however, would not alter the results 
significantly.) 

We now seek to relate the steady disc accretion rate, M, to the 
observed accretion rate on to the stellar surface. Mace - This depends 
on the flow of material across the gap i nduced by the planet , and is 
not yet fully understood. The models o f iLubow & d' Angeld |2003) 
suggest that the accretion rate past the planet (into the inner disc) is 
typically 10-25% of the accretion rate on to the planet. The accre- 
tion rate on to the planet, Mp, in turn depends on the planet mass 
Mp. This interaction is still the subject of some debate , but for sim- 
plicity we adopt the approximate scaling derived bv iLubow et al.l 
i l 19991) 
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and, following iLubow & d' Angeld booel) , set 
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Consequently, the observed accretion rate and the disc mass are 
related by 



Mace ^ 0.25 
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which is a function only of the planet mass, viscosity and disc size. 
For simplicity, we fix R„^t = lOOAU and consider a range of planet 
masses and viscosities. We consider planet masses ranging from 
IMjup (approximately the minimum mass required to open a gap 
in the disc) to 12M]„p (approximately the deuterium-burning limit), 
and viscosity parameters ranging from a = 0.01-0.001. 

Fig[5]shows the range of parameters allowed by this model, as 
well as that of the photoevaporation model. As discussed in Section 
14.11 objects appear as inner hole sources in the photoevaporation 
model only when the accretion rate falls below the wind rate, and 
unless the disc is unusually inviscid, or unusually large, this implies 
a maximum disc mass of 0.005Mo. In principle it is possible to 
observe a "photoevaporated" inner hole with a detectable accretion 
rate (as in the fourth snapshot of Fig[3]l, but the window of oppor- 
tunity for such a detection is very short and the majority of such 



Properties of inner hole sources 




GoKu Tou/f" Passive 



Md / Mo 

Figure 5. Disc mass versus accretion rate for different models of inner hole 
sources. The solid lines show the result of the "planef model: solid lines 
show the results for a = 0.01; dashed lines denote a = 0.001. In each case 
the upper and lower lines show the results for the minimum (IMjup) and 
maximum (12Mj„p) planet masses respectively. The hatched region there- 
fore represents the range of parameters allowed by the planet model. The 
solid grey region shows the range allowed by the photoevaporation model 
but, as noted in the text, most such inner hole sources are not expected to 
have detectable accretion. Also shown are data for several known inner hole 
sources, as described in the text: circles represent detections, while triangles 
and arrows denote upper limits. The error bar in the upper left corner shows 
the typical systematic error in the observational derivation of these parame- 
ters. 



objects have no significant accretion. By contrast, planet-induced 
inner holes tend to show significant accretion rates, and can occur 
over a range of disc masses. For a given disc mass a more massive 
planet results in a lower accretion rate, and a lower viscosity results 
in a similarly lower accretion rate. 

Also shown in Fig[5] are the properties of several well- 
studied transition objects: GM Aur, DM Tau, TW H y a and 
CoKu Tau/4. These data a re tak en from ICalvet et alj j2002l 
TW Hya) and Calyetetal] ilOOj. GM Aur & DM Tau), with 
disc m asses for sources in Taurus -Auriga taken from the sur- 
vey of [Andrews & Williams! ( |2005|) . CoKu Tau/4 is a spectro- 
scop ic WTT with no ev i dence of on-going accretion so, follow- 
ing S icilia-Aguilar et al] ( l2006h . we assign an upper limit to the 
accretion rate of lO^'^Mgyr"'. In addition, we have included the 
locus of the "passive discs" identified in the survey of binary 
objects conducted by iMcCabe et al.l ( |2006|) . This survey identi- 
fied six inner hole sources, four of which are in Taurus- Auriga. 
[Andrews & Williamll l l2005b derived disc masses of ^ 0.0005- 
O.OOIMq for these objects, but were not able to resolve the indi- 
vidual components of these binary systems. Consequently we treat 
these measurements as upper limits to the disc masses, as they 
essentially measure the total disc masses of the binary pairs. All 
of these objects are WTTs, so we again assign an upper limit of 
lO^'-Moyr-'. 

This toy model obviously has uncertainties of factors of a few, 
but despite this Fig[5] suggests that observations of disc masses and 
accretion rates potentially provide a powerful tool for discriminat- 
ing between different models. GM Aur, DM Tau and TW Hya fall 
in the region predicted by the planet model, w hile CoKu Tau/4 and 
the objects identified bv lMcCabe et al.l ( l2006h are more readily ex- 
plained by photoevaporation. The reason that this simple method 
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can discriminate so strongly between tiiese models is easily un- 
derstood. We select objects based on the "transitional" appearance 
of their SEDs only, but then consider parameters that relate to the 
global disc structure. An embedded planet affects the SED dramat- 
ically, but has only a modest (factor of ~ 10-100) effect on the ac- 
cretion rate and little effect on the total disc mass (as the most of 
the disc mass resides at large radii). By contrast, photoevaporative 
clearing requires significant changes in the global disc properties, 
and therefore occupies a markedly different region of parameter 
space in Fig[5] 

We note, however, that our toy planet model does not con- 
sider the formation of the planet, but merely assumes its presence. 
More detailed analysis of this problem can provide additional in- 
sight, and consideration should also be given to the evolution of the 
planet-disc system over ~Myr time-scales. Additional parameters, 
such as stellar ages and masses, can also be used to constrain mod- 
els further and more detailed observations of disc properties can 
add valuable additional constraints (e.g. lEisner et al.ll2006) . How- 
ever, our results suggest that even simple, easily observed parame- 
ters, such as disc masses and accretion rates, can provide a strong 
discriminant between different models of disc clearing. 



6 SUMMARY 

In this paper we have studied the dynamics of dust and gas during 
protoplanetary disc clearing. We have constructed two-fluid numer- 
ical models, in which the gas is subject to viscous accretion and 
photoevaporation and the dust undergoes both radial diffusion and 
advection. We find that small grains (< 10pm) grains remain well- 
coupled to the gas throughout, while larger grains undergo signifi- 
cant radial migration over the duration of the clearing phase. This 
suggests that mm grains in the outer disc must be replenished on 
time-scales of 10''-10'yr, presumably through coUisional growth. 
Once the photoevaporative wind opens a gap in the gas disc the 
dust is rapidly removed from the inner disc. As the outer gas disc is 
cleared by the wind the grains are "swept up" by the steep pressure 
gradient at the inner edge of the gas disc, resulting in a dust-rich 
ring near to the inner disc edge. This may result in increased rates 
of grain collision and growth, with the potential to form planetesi- 
mals and/or a debris disc. We have described in detail the observa- 
tional appearance of discs during the clearing phase in our models, 
and find that several observed discs agree well with the predictions 
of our model. Lastly, by considering a simple model of a planet em- 
bedded in a disc we have outlined how observations of the masses 
and accretion rates of transitional discs can provide a means of dis- 
criminating between different models of disc clearing. 
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APPENDIX A: WIND PROFILES 

The radial mass-loss profiles of the diffuse- and direct-field photo- 
evaporative winds are defined as follows. In the "diffuse" regime, 
where the inner disc is optically thick to ionizing photons, the 
wind is driven by recombinations in the disc atmosphere . We adopt 
the w ind profile from the hydrodynamic simulations of iFont et al.l 
( 12004 kindly provided in numerical form by Ian McCarthy): 



SdiffiiscC^) = 2no(R)ui{R)ijmii . 



(Al) 



Here no(^) is the density at the base of the ionized layer, U](R) the 
wind launch velocity and /j the mean molecular weight of the ion- 
ized gas (taken to be yu = 1.35 throughout). The base density profile 
is given by the fitting form 



no(R) 



I , ^25/2 

(f) ^(f) 
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where the base density at R„ = GM,/cl (where 
the sound speed of the ionized gas) is given by 



(A2) 



10km s is 
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[4naBRl, 

Here <I> is the stellar ionizing flux (in units of photons s"') and qtb is 
the Case B recombination coefficient for atomic h ydrogen at lO^'K, 
which has a value of as = 2.6 x lO-'^'cm^'s"' jcoxll2000 f). and 
the constant Ci - 0.14. The numerical simulations show that the 
launch velocity profile is well-matched by the numerical form 



U](R) = CjAexp 



, R 



0.1 



A -0.1 

R., 



R^O.IR, (A4) 



where the numerical constants have the values A = 0.3423, B = 
-0.3612 and D = 0.2457. For radii smaller than Q.IR^ the wind 
rate is set to zero. 

In the "direct" regime, where the inner disc is optically thin to 
ionizing photons, the wind is instead driven by direct illumination 
of the inner disc edge. In this case we ado pt the wind profile from 
the simulations of I Alexander et al.l ( l2006al) . The base density at the 
inner disc edge is set by ionization balance (essentially a Stromgren 
criterion), and the hydrodynamic simulations show that the follow- 
ing power-law form provides a good fit to the wind profile: 



\Ana^{HIR)R?Jt) 
R 



1/2 



(A5) 



R,At) 



R>R„ 



Here /?in(0 is radius of the inner disc edge, which increases as the 
wind disperses the outer disc. The fitting constants cJ3 and a vary 



^ Note that the consta nt defined as C2 here is denoted by "CD" in 
[Alexander et afj i2006d) . 



14 R.D. Alexander & P.J. Armitage 

somewhat with the disc aspect ratio HjR, but for simphcity we 
adopt the values derived for H/R = 0.05: C = 0.235 and a = 2.42. 



